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ABSTRACT 

We present X-ray proper-motion measurements of the forward shock and reverse-shocked ejecta 
in Tycho's supernova remnant, based on three sets of archival Chandra data taken in 2000, 2003, 
and 2007. We find that the proper motion of the edge of the remnant (i.e., the forward shock 
and protruding ejecta knots) varies from 0".20yr _1 (expansion index m = 0.33, where R = t m ) to 
0".40yr~ 1 (m = 0.65) with azimuthal angle in 2000-2007 measurements, and 0".14yr -1 (m = 0.26) 
to 0".40yr~ 1 (m = 0.65) in 2003-2007 measurements. The azimuthal variation of the proper motion 
and the average expansion index of ~0.5 are consistent with those derived from radio observations. 
We also find proper motion and expansion index of the reverse-shocked ejecta to be 0". 21-0". 31 yr _1 
and 0.43-0.64, respectively. From a comparison of the measured m-value with Type la supernova 
evolutionary models, we find a pre-shock ambient density around the remnant of < 0.2 cm -3 . 

Subject headings: ISM: individual object (Tycho's SNR) — ISM: supernova remnants — shock waves 
- X-rays: ISM 



1. INTRODUCTION 

A supernova (SN) event which Danish astronomer Ty- 
cho Brahe discovered in 1572 is now visible at many wave- 
lengths as Tycho's SN remnant (SNR). The SN event has 
been classified as a normal Type la (i.e., a thermonuclear 
explosion of a white dwarf star in a close binary) based 
on the shape of its light curve (e.g., Baade 1945; Ruiz- 
Lapuente 2004) and the X-ray properties (Badenes et 
al. 2006). This is recently confirmed from its light-echo 
spectrum, which reached the earth about 436 yr after 
Tycho Brahe observed the direct light from the SN event 
(Krause et al. 2008). Given the fact that Type la SNe 
have been used to determine cosmological parameters 
such as the Hubble constant and the Dark Energy equa- 
tion of state (e.g., Hamuy et al. 1996; Riess et al. 1998), it 
is fundamentally important to understand the explosion 
mechanism and the progenitor system. Tycho's SNR is 
a very suitable target for such a study. First, its ejecta- 
dominated X-ray emission allows us to constrain the ex- 
plosion mechanism models as was recently demonstrated 
by Badenes et al. (2006), who concluded that the explo- 
sion producing Tycho's SNR was most likely a delayed 
detonation; also, the surviving binary companion from 
Tycho's SNR may have been identified (Ruiz-Lapuente 
et al. 2004; Ihara et al. 2007), though the association 
with the SN event is still controversial (Kerzendorf et al. 
2009). The chemical abundances of the candidate com- 
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panion were recently measured to be consistent with the 
Galactic trend except for an overabundant Ni/Fe ratio 
which suggests contamination by SN ejecta (Hernandez 
et al. 2009). Recent detection of Cr and Mn Ka lines in 
the X-ray spectrum of the SNR (Tamagawa et al. 2009; 
Yang et al. 2009) also provides an opportunity to con- 
strain the metallicity of the progenitor star (Badenes, 
Bravo, & Hughes 2008). 

Another important aspect of Tycho's SNR is the pres- 
ence of non-thermal X-ray synchrotron emission pro- 
duced by accelerated particles (e.g., Pravdo & Smith 
1979). Recent Chandra images have revealed that the 
non-thermal emission is confined to a very thin shell just 
behind the forward shock (Hwang et al. 2002; Bamba 
et al. 2005; Warren et al. 2005). The thin filamentary 
structure of non-thermal emission is now known to be a 
common feature in many young SNRs, while its nature is 
the subject of considerable theoretical work on particle 
acceleration (e.g., Pohl et al. 2005; Ballet 2006; Cassam- 
Chena'i et al. 2007). Furthermore, this SNR is a possible 
TeV gamma-ray source, although a firm detection of TeV 
gamma-ray emission is not yet available (e.g., Volk et al. 
2008). 

Tycho's SNR potentially provides us with information 
on SNR evolution, since the degree of deceleration of the 
forward shock can be inferred by combining the current 
expansion rate with the accurately known age. Expan- 
sion measurements have been performed at optical, radio, 
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cused on the northwestern (JN W), northern (N), and east- 
ern (E) Ho filaments. The expansion index (m; R oc t m , 
where R and t are the radius and the age of the remnant) 
for these filaments was measured to be 0.32-0.41 (Kam- 
per & van den Bergh 1978), indicating significant decel- 
eration from free expansion (m = 1). The radio obser- 
vations revealed that the expansion index varies between 
0.2 and 0.8 with azimuthal angle (e.g., Reynoso et al. 
1997). The slow expansion (m ~ 0.2) seen at the E rim 
suggests an interaction between the forward shock and a 
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dense cold cloud. Such an interaction was confirmed by 
the subsequent discovery of a dense cloud toward the E 
rim (Reynoso et al. 1999; Lee et al. 2004). The mean m 
value was calculated to be ~0.47, suggesting that over- 
all, this SNR is still evolving toward the Sedov phase 
(m = 0.4). Although the mean value of m for radio ob- 
servations is higher than that for optical observations, 
optical and radio measurements of individual features 
are consistent with each other. In X-rays, the only X- 
ray expansion measurement to date was performed using 
ROSAT data with moderate angular resolution of ^8" 
half-power diameter (Hughes 2000). Interior features of 
the remnant show consistent expansion indices between 
X-ray and radio observations. In contrast, the expan- 
sion index for the forward shock measured in X-rays, 
m ~ 0.71, is higher than that for radio. Since the radio 
and X-ray images are well correlated with each other, the 
expansion rate discrepancy has been puzzling. 

Here, we present X-ray proper-motion measurements 
based on archival Chandra data with superior angular 
resolution of ^0.5" half-power diameter taken in three 
epochs between 2000 and 2007. These Chandra data al- 
low us to measure proper motions of the SNR edges as 
well as the reverse-shocked ejecta inside the remnant. 

2. OBSERVATIONS AND DATA REDUCTION 

Chandra observed Tycho's SNR in three epochs as 
summarized in table [TJ The third epoch comprises two 
observations taken on April 23 and 26 in 2007 (ObsID 
7639 and 8551; PI: Hughes). The exposure time of Ob- 
sID 8551 was significantly shorter than that of ObsID 
7639. For simplicity, we use only the longer observation 
(ObsID 7639) in our analysis. The information of the 
Chandra observations used in this paper is listed in. The 
time differences between the first and the third observa- 
tions (i.e., 2000-2007) and between the second and the 
third (i.e., 2003-2007) are 6.60 yr and 3.99 yr, respec- 
tively. The observations in 2003 and 2007 were specif- 
ically intended to allow proper-motion measurements: 
both have the same pointing direction, roll angle, and 
exposure time. The Advanced CCD Imaging Spectrom- 
eter (ACIS)-I was used for these two observations. On 
the other hand, the 2000 observation was not intended 
for this purpose, and utilized the back-illuminated ACIS- 
S3 detector. However, the 2000 observation is still use- 
ful to measure proper motions. Therefore, we use the 
2000 observation, which in fact provides slightly tighter 
constraints on proper motion due to the longer inter- 
val between observations. We reprocessed all the level- 1 
event data, applying the standard data reduction^ with 
CALDB ver. 3.4.0. The net exposure times for the 2000, 
2003, and 2007 observations are 48.9 ks, 145. 7ks, and 
108.9 ks, respectively. In order to make use of the best 
position accuracy of the ACIS, we removed the effects of 
pixel randomization. 

Figure [T] shows a three-color Chandra image obtained 
from the 2003 observation. Red corresponds to low- 
energy (0.3-1.0keV, mostly Fe L), green to middle- 
energy (1.0-4.0 keV, mostly Si+S+Ar+Ca K), and blue 
to high-energy (4.0-8.0 keV, mostly Fe K blend plus the 
continuum) bands. The remnant is outlined by thin 
blue (nonthcrmally-dominated) filaments, while it con- 
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tains substantial green (thermally-dominated) fragmen- 
tation inside the blue filaments. The blue filaments are 
identified as the forward shock (e.g., Warren et al. 2005), 
while the fragments are the reverse-shocked ejecta (e.g., 
Velazquez et al. 1998). We measure proper motion of 
both the forward shock and the shocked ejecta, by com- 
paring vignetting-corrected 1.0-8.0 keV images obtained 
in the three epoches. 

3. PROPER MOTION MEASUREMENTS 

Before measuring proper motions, we try to register 
the three images by aligning the positions of several ob- 
vious X-ray point sources surrounding the remnant. We 
determine their positions in each epoch, by employing 
wavdetect software included in CIAO ver. 4.0. We de- 
tect nine sources around the SNR in the 2003 and 2007 
images, while only five are detected along the NW pe- 
riphery in the 2000 image. The smaller number of sources 
in the 2000 image is partly due to the fact that the south- 
ern portion of the SNR was not covered by ACIS-S3 chip, 
and partly due to the shorter exposure time. If we align 
only the sources concentrated in the small NW area, it 
may cause misalignment for the rest of the image. There- 
fore, we do not correct the registration between the 2000 
and 2007 observations, and instead incorporate into our 
proper motion measurement for 2000-2007 a systematic 
uncertainty of 0".4, which is the 1-sigma absolute astro- 
metric accuracy for the ACIS camera reported by the 
Chandra calibration tearr0. This systematic uncertainty 
seems to be conservative, because the rms residual cal- 
culated for the five point sources detected in both im- 
ages is only 0".13. On the other hand, we use the point 
sources to register images between 2003 and 2007. Since 
no optical counterpart was identified within a 1" radius 
for any of the point sources based on any catalogs con- 
tained in the VizieR catalog servic^J we assume that 
none have undergone proper motion between 2003 and 
2007. We aligned the coordinates of the images so that 
the point sources have the same positions, by using ccmap 
and ecsetwes in IRAF IMCOORDS. In the alignment 
procedure, four parameters for the second-epoch image 
are allowed to vary: x and y shifts, rotation, and pixel 
scale. The best-fit parameters show differences for both 
x and y shifts of less than 0".2, no modifications of pixel 
scales, and rotation of 3', respectively. After the align- 
ment, the rms residual for the positions of the nine point 
sources between 2003 and 2007 images is ^0".29, which 
will be taken as the systematic uncertainty in our proper- 
motion measurements for 2003-2007. 

Figure [2] left and right show the difference images in 
the 1.0-8.0keV energy band for 2000-2007 and 2003- 
2007, respectively. The subtraction is performed after 
each image was registered and smoothed by a Gaussian 
kernel of a = 0".492 (i.e., 1 ACIS pixel) which is ap- 
proximately the angular resolution of Chandra's X-ray 
telescope. There are some artifacts in the difference im- 
ages: the southern dark (negative) region seen in the 
2000-2007 difference image which is due to the fact that 
the 2000 image has no exposure there, the structures 
running from NW to SE and from NE to SW caused by 
the ACIS chip gaps, and stripes running from NE to SW 
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caused by bad columns. Apart from these artifacts, we 
can clearly see motion of the X-ray features as adjacent 
black and white features for the forward shock (the edge 
of the remnant) and the reverse-shocked ejecta (knotty 
features inside the remnant). We quantify the proper 
motion of these features in the following sections. 

3.1. SNR Edges 
3.1.1. Forward Shock 

In order to measure proper motions accurately for the 
forward shock, we generate one-dimensional X-ray pro- 
files along the apparent direction of shock motions, fol- 
lowing Katsuda, Tsunemi, & Mori (2008). These one- 
dimensional X-ray profiles are extracted from the 39 
small boxes (25"x25") indicated in Fig.Q] The position 
angle of each box is chosen such that it is as perpendic- 
ular as possible to the shock front (by eye). The choice 
of position angle does not affect the proper motion mea- 
surement within at least ±6 degrees difference, which we 
checked for several boxes. In order to include the cali- 
bration uncertainty of the ACIS effective area, we add a 
f03 systematic error term in quadrature to the statistical 
error for each data point. Each profile is binned by 0".5 
along the radial direction of the SNR. An example pro- 
file is shown in Fig. [3l in which the shift due the motion 
of the shock is readily apparent. Shifting the first-epoch 
profile and calculating the x 2 values from the difference 
between the two profiles at each shift position, we search 
for the best-matched shift position for each box. Fig- 
ure |4] shows the x 2 distribution as a function of shift 
position for the example profiles in Fig. |3] By applying 
a quadratic function for the x 2 profile, we measure the 
best-shift position where the minimum x 2 value occurs. 
In this way, we obtain two (2000-2007 and 2003-2007) re- 
sults for each box. The reduced-x 2 value, which is nearly 
equal to x 2 value divided by the number of bins in each 
box (here 36), falls in the ranges 0.9-10.4 for 2000-2007 
and 0.9-2.6 for 2003-2007 at the best-matched shift posi- 
tions. Two neighboring boxes (azimuth angles = 78° and 
85° counterclockwise from N) for the 2000-2007 measure- 
ments have large (>3.7) reduced-x 2 values. This is partly 
due to the fact that these two include the brightest parts 
among the entire SNR edge, resulting in smaller statis- 
tical errors than those of the others, and partly due to 
the fact that the angular resolution around the two boxes 
are significantly different between 2000 and 2007 obser- 
vations. At the location of the two boxes, we simulate 
point spread functions (PSFs) for 2000 and 2007 images 
by using chartP^I and mar J 1 ^! software, and find that the 
shape of the PSF for the 2000 image is elongated in the 
direction of shock motion, while that for the 2007 image 
is not. If we fit the PSFs with a Gaussian kernel, we 
obtain a ~ 1".7 and a ~ 0".9 for the 2000 and 2007 im- 
ages, respectively. Therefore, we smooth the 2007 image 
by a Gaussian kernel of a — 1".476 (i.e., 3 ACIS pixels), 
such that it has the same image resolution as that of the 
2000 image. In this way, we obtained reasonable fits for 
all the box regions (reduced-x 2 < 2.6). 

Also notable in Fig. [3] is possible flux variations among 
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the three epochs. We find that X-ray flux variations 
among the three epochs are within ~10% for all the 
boxes. The 10% flux changes are larger than the ACIS 
effective area calibration uncertainty (5%). Consider- 
ing that the magnetic field strength in Tycho's SNR is 
thought to be in the milligauss range (e.g., 0.1-10 mG, 
Reynolds k Ellison 1992; 0.03-0.4 mG, Warren et al. 
2005), year-scale flux variations of the synchrotron emis- 
sion seen in RXJ11713. 7-3946 (Uchiyama et al. 2007) 
or Cas A SNRs (Uchiyama & Aharonian 2008) may be 
expected in Tycho's SNR as well. However, the possi- 
ble flux variations in Tycho's SNR are small compared 
with those seen in these other SNRs, where some knots 
showed more than 100% flux variations. Further detailed 
investigations are needed to determine whether the flux 
variation is due to the flux change of the synchrotron 
emission, and is beyond the scope of this paper. 

3.1.2. Ejecta Knots Protruding beyond the Forward Shock 

There are breakout regions along the periphery of the 
SNR, as indicated by the small annular segments with 
triangular marks in Fig. [T] In these regions, the SNR 
edges show fragmented features which look like Rayleigh- 
Taylor fingers. It is already shown that this fragmenta- 
tion is formed by ejecta knots protruding beyond the for- 
ward shock, based on spectral analysis of XMM-Newton 
data (Decourchelle et al. 2001) as well as principal com- 
ponent analysis of Chandra data (Warren et al. 2005). 
Since these ejecta knots show completely different mor- 
phology from the filamentary structure of the forward 
shock, we cannot define box regions relative to the shock 
front. We thus use two-dimensional images instead of 
one-dimensional profiles to measure the proper motion. 
We follow the method employed by DeLaney & Rud- 
nick (2003), who measured proper motion of the forward 
shock in Cas A. We first generate polar-coordinate im- 
ages for the three observations. Then, radially shifting 
the polar-coordinate image from 2000 or 2003 by steps 
of 0".5 with respect to 2007 image, we search for the 
best-matched shift positions. To quantitatively obtain 
the best-shift position, we calculated the C value at each 
shift position, where C is expressed as (Cash 1979; Vink 
2008), 

C = — 2\nP = — 2 y ^(njnmj — m, — lnnj), 

i 

with rii being the observed counts in a pixel i and be- 
ing the model counts in the same pixel. We here assume 
that the 2007 observation is the model image. Using the 
C value rather than x 2 is necessary when the number of 
counts is so low that the shape of the Poisson distribution 
becomes significantly different from that of the Gaussian 
distribution. In fact, this is the case for a number of 
pixels in the data used here. Since the analyzed areas 
(i.e., the small annular segments in the Fig. [T]) include 
regions outside the remnant, they contain a significant 
number of pixels with zero counts. Although we cannot 
sum pixels with rii (or Wj) = according to the above 
equation, we have to sum the same numbers of pixels 
for each shift position. To this end, we assume the ob- 
served and model counts in the pixels having zero counts 
to be unity. We check the validity of this assumption, 
by applying the two-dimensional method for several for- 
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ward shock regions and confirming that the proper mo- 
tion value obtained from the one- and two-dimensional 
methods are consistent with each other. As an exam- 
ple, the C values calculated for the yellow box in Fig. [1] 
is shown as a function of shift position in Fig. |4] (with 
dashed lines), from which we can see that the best-shift 
position from the two-dimensional method (C value) is 
the same as that from the one-dimensional method (x 2 
value). The best-shift position for each area is estimated 
by fitting a quadratic function for each profile. In esti- 
mating uncertainties, AC (= C m - m — C) can be treated 
as Ax 2 (Cash 1979). In this way, we obtained two sets 
(2000-2007 and 2003-2007) of proper motion measure- 
ments for each small annular segment. 

3.1.3. Results for the SNR Edge 

Table [5] gives the proper motion, expansion rate, and 
expansion index (which are derived from "the expansion 
rate" x "the age of the remnant, i.e., 430 or 432 yr for 
2000-2007 or 2003-2007, respectively") for all areas in- 
cluding the forward shock (small boxes in Fig. [IJ and 
ejecta knots (small annular pieces in Fig.[T]). We also list 
the mean values which are simply the sum for all areas 
(both 2000-2007 and 2003-2007 measurements) divided 
by the number of areas. To calculate the expansion rate 
and expansion index, we need to determine the expan- 
sion center. We take the geometric center (J2000) of 
RA=00 h 25 m 19 s .4, DEC = +64°08'13".98 determined by 
minimizing the ellipticity of the forward shock (Warren 
et al. 2005). Figure [5] plots the values as a function of 
azimuthal angle. Black and red points represent 2000- 
2007 and 2003-2007 measurements, respectively. Data 
points with filled boxes are responsible for the forward 
shock and open circles for the ejecta knots. The errors 
quoted in Fig. [5] represent statistical 1-sigma uncertain- 
ties. The magnitude of the systematic uncertainties are 
indicated in the figure as black and red bars, and are 
generally ~3 times larger than the statistical uncertain- 
ties. If both statistical and systematic uncertainties are 
considered, the two results for each area are all consis- 
tent with each other except for the E rim (70°-90° in 
azimuthal angle) where we see indications of significant 
shock deceleration. The variation of the proper motion 
ranges from 0".14yr _1 at the E rim to 0".40yr _1 at the 
western (W) rim in 2003-2007 measurements, but is rela- 
tively more modest variation at 0".20yr _1 to 0".40yr _1 
in 2000-2007 measurements. 

3.2. Reverse-Shocked Ejecta 

We now describe our measurement of the proper mo- 
tion of the reverse-shocked ejecta. Warren et al. (2005) 
located the position of the reverse shock at ~3' radius, 
which is the peak of the radial profile of the Fe K shell 
line. They also reported that the contact discontinu- 
ity between the reverse-shocked ejecta and the swept-up 
interstellar medium (ISM) is located around 3'. 67-4'. 5 
radius. We thus focus on an annular region between the 
reverse shock (at 3' radius) and the contact discontinuity 
(at 3'. 71-4' radius depending on azimuthal angle, based 
on Fig. 4 in Warren et al. 2005). To investigate azimuthal 
variations, we divide the annular region into five sectors 
as shown in Fig. [1] The five sectors represent NE, south- 
east (SE), southwest (SW), W, and NW defined as: NE 



(0°-60° in azimuthal angle), SE (132°-180°), SW (180°- 
240°), W (240°-300°), NW (300°-360°). Note that we 
exclude the E sector from our analysis. In these regions, 
Warren et al. (2005) could not identify the reverse shock 
due to the very weak Fc K line emission there, nor do we 
see clear evidence of ejecta motion. 

As mentioned in Sec. 2, the emission from the reverse- 
shocked ejecta appears as fragmented structures and 
is completely different from the filamentary structures 
of the forward shock. Therefore, we used the two- 
dimensional method used for measuring the proper mo- 
tion of the ejecta knots (see Sec. 3.1.2). The C-value cal- 
culated for the example sector shown in yellow in Fig.[T]is 
plotted as a function of shift position in Fig. [6] Table [3] 
lists the proper motion, expansion rate, and expansion 
index for individual regions with their statistical uncer- 
tainties as well as their mean values. Figure [7] shows 
the results as a function of azimuthal angle. Black and 
red points represent 2000-2007 and 2003-2007 measure- 
ments, respectively. Errors with the best-fit values rep- 
resent only statistical uncertainties, while the systematic 
errors are shown as black and red bars. The azimuthal 
variations for the proper motion range from 0".21 yr -1 at 
the NW sector to 0".31yr" 1 at the SW sector, although 
they do not show significant variation within their com- 
bined statistical and systematic uncertainties. 

It is of interest to investigate the X-ray energy depen- 
dency of the expansion, especially for the reverse-shocked 
ejecta whose emission lines show distinct radial distribu- 
tions: the radial profiles of Si K and S K line emission 
peak outside that of the Fe K line emission (Hwang & 
Gotthelf 1997; Decourchelle et al. 2001). Furthermore, 
Dopplcr velocities of the Si K and S K lines were recently 
reported to be larger than that of Fe K line (Furuzawa et 
al. 2009; Hayato et al. in preparation), indicating higher 
velocity of the Si and S ejecta than of Fe. Therefore, we 
performed expansion measurements of the shocked ejecta 
for selected energy bands. Unfortunately, we could not 
derive meaningful results for a narrow energy band con- 
taining the Fe K line due to insufficient statistics. On 
the other hand, we find that results for the energy bands 
of 0.75-1.43 keV (dominated by the Fe L lines) and 1.63- 
4.1 keV (dominated by Si and S K lines), which are not 
shown here, are both consistent with those derived from 
the 1.0-8.0keV band, as we expect given that images in 
these energy bands are very similar (Hwang & Gotthelf 
1997; Decourchelle et al. 2001). 

4. DISCUSSION 

4.1. Comparison among Radio, Optical, and X-Ray 
Measurements 

The proper motion along the edge of the SNR (i.e., for- 
ward shock and ejecta knots) varies with azimuthal an- 
gle. In general, the variation is quantitatively consistent 
with radio measurements in which slower motions are 
observed at the E rim and faster motions at the W rim 
(0".17yr- 1 -0".31yr~ 1 by Strom et al. 1982; 0".14yr~ 1 - 
0".42yr _1 by Reynoso et al. 1997), Our results are also 
consistent with the optical measurements of the NW, N, 
and E rims (Kamper & van den Bergh 1978). The aver- 
aged expansion index for all areas is ~0.52. This value is 
consistent with or close to those determined by radio ob- 
servations: 0.47±0.07 (Strom et al. 1982), 0.462±0.024 
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(Tan & Gull 1985), 0.471±0.028 (Reynoso et al. 1997). 

4.2. Interpretation of the Azimuthal Variation 

The shock motion in the E rim is the slowest along 
the entire rim. This agrees with the radio proper mo- 
tion measurements and further supports the idea that 
the forward shock is interacting with a dense cold cloud 
in this region (Strom 1983; Reynoso et al. 1999; Lee et 
al. 2004). The interaction is suggested to be very recent 
(within the last 50 yr) and strong (the cloud density of 
160-325 cm -3 ) from radio observations (Reynoso et al. 
1999). The X-ray data show a hint of deceleration along 
eastern rim (especially at 70°-90° in azimuthal angle). 
This possible deceleration, if true, would strongly sup- 
port the view that the Tycho's SNR has been recently 
interacting with the dense cloud along the E rim. 

In Fig. [5] left, we also see a possible variation of the 
proper motion of the forward shock along the W-NW 
rim. The value of proper motion is generally propor- 
tional to the shock radius in this region (see, Fig. [5] or 
Fig. 4 in Warren et al. 2005), resulting in fairly a con- 
stant expansion rate throughout the region as shown in 
Fig. [3] right. The gradual increase of the proper motion 
from the NW rim to the W rim suggests that the ambi- 
ent density along the NW rim is denser than that along 
the W rim. Since there is no evidence in the radio of a 
dense cloud around the NW rim (Reynoso et al. 1999; 
Lee et al. 2004) , a recent strong shock deceleration does 
not seem to be the case. Therefore, we suggest that the 
forward shock has been constantly decelerated in these 
regions. The slower shock motion as well as the smaller 
shock radius in the NW direction most likely reflect a 
slightly denser ambient medium in this direction than 
that in the W direction. We can safely assume pressure 
equilibrium behind the forward shock (i.e., nowf, where 
no is the ambient density and v s is the shock velocity) 
along the NW rim to the W rim. Using the proper mo- 
tions of ~0".26yr~ 1 at the NW rim and -0".34yr _1 at 
the W rim, we find that the ambient density ratio of NW 
rim to W rim is *~2. The higher ambient density toward 
the NW direction should yield a stronger (i.e., slower) re- 
verse shock in this region. This is indeed consistent with 
the brighter X-ray intensity of the ejecta in the NW sec- 
tor than in the W sector (see Fig. [T|) and the possibility 
that the reverse-shocked ejecta may be moving slower in 
the NW sector than in the W sector (though this is not 
significant within the systematic uncertainties). 

4.3. Ambient Density 

We estimate the ambient density around Tycho's SNR, 
based on a comparison of the measured expansion index 
with theoretical SNR evolutionary models. Dwarkadas 
& Chevalier (1998) used one-dimensional simulations to 
investigate SNR evolution for power-law and exponen- 
tial ejecta density profiles. They predicted that the ex- 
ponential profile shows a density curve increasing from 
the reverse shock to the contact discontinuity, while 
the power-law profile shows a decrease. The density 
profile expected for the exponential profile qualitatively 
matches that observed in Tycho's SNR (Hwang & Got- 
thelf 1997). Furthermore, the exponential profile is 
a better approximation than the power-law profile for 
the SN ejecta in detailed explosion models. Dwarkadas 



& Chevalier (1998) thus concluded that the exponen- 
tial profile is more suitable for Tycho's SNR than the 
power-law profile. Dwarkadas (2000) extended the evo- 
lutionary model with the exponential ejecta density pro- 
file to two dimensions, including the effect of Rayleigh- 
Taylor fingers acting at the contact discontinuity. We 
compare our measured expansion index for the forward 
shock (i.e., SNR edge) with that expected from the 
two-dimensional model examined in a uniform ambient 
medium. In this model, the evolution of the expansion in- 
dex of the forward shock can be scaled by a characteristic 

time, 248E~ 1 1/2 (M ci /M ch ) 5 / 6 n~ 1/3 yr (e.g., Dwarkadas 
& Chevalier 1998), where E$i is the explosion energy in 
units of 10 51 ergs, Mch is defined to be 1.4 M Q , and no 
is the ambient H density in units of cm -3 . We exclude 
the E rim where the cloud-shock interaction is obvious 
(60°-100° in azimuthal angle), since the model is calcu- 
lated assuming spherical symmetry. Thus, we use an av- 
erage m-value for the SNR edge of ~0.54. This m-value 
requires an age equal to the characteristic time accord- 
ing to Fig. 1 in Dwarkadas (2000). For the age of Ty- 
cho's SNR (432 yr), we derive the ambient density to be 

-0.2£; 5 " 1 3/2 (A/ ej /A/ C h) 5/2 cm- 3 . This density is at least 
three times lower than that required from a hydrody- 
namic model to reproduce structures of the ionization age 
and the electron temperature seen in the X-ray line emis- 
sion (Dwarkadas & Chevalier 1998). Also, the density 
is about five times lower than that preferred by detailed 
comparison between hydrodynamical Type la SN models 
and observed line emission in X-ray spectra (Badenes et 
al. 2006). On the other hand, this density is well within 
the upper limit for the ambient density of 0.3 cm -3 (with 
an uncertainty of factor ~ 2) based on the contribution of 
the thermal emission from the shocked ambient medium 
in the Chandra spectra (Cassam-Chenai et al. 2007). It 
is also consistent with both the previous optical measure- 
ment of 0.3 cm -3 for the NE filament (Kirshner, Winkler, 
& Chevalier 1987) and the upper limit of 0.4 cm -3 based 
on the gamma-ray flux (Volk et al. 2008). Therefore, 
comparisons of hydrodynamical models with the X-ray 
line emission agree in the need for higher densities than 
the other estimates. As Dwarkadas & Chevalier (1998) 
noted, it might be worth considering a more complex cir- 
cumstellar and/or ejecta structure in the hydrodynami- 
cal model to explain the X-ray line emission. 

We note that the effects of efficient particle acceler- 
ation are not considered by the models above, but are 
now believed to influence the forward shock in this SNR 
(Warren et al. 2005; Cassam-Chenai et al. 2007). Effi- 
cient particle acceleration would yield high compression 
ratios behind the shock well in excess of 4, resulting in 
stronger shock deceleration than if efficient particle ac- 
celeration were absent. 

There are a number of authors who considered parti- 
cle acceleration effects on the SNR evolution (e.g., De- 
courchelle, Ellison, & Ballet 2000; Ellison et al. 2007). 
Ellison et al. (2007) showed in Fig. 2 in their paper that 
the model with e re i i2 63% roughly matches the radius 
ratios among the forward and reverse shocks, and the 
contact discontinuity observed in Tycho's SNR, where 
e rc i is defined as the percentage of energy flux crossing 
the shock that ends up in relativistic particles. This 
model predicts that the radius and velocity of the forward 
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TABLE 1 

Chandra OBSERVATIONS USED IN THIS PAPER 



ObsID 


ObsDate 


Array 


Coordinate (J2000) 




Exposure Time (ksec) 


115 


2000-9-20 


ACIS-S 


00 h 25 m 07 s .34, 64°09'42" 


7 


48.9 


3837 


2003-4-29 


ACIS-I 


00 h 25 m 21 s .48, 64°08'13" 


9 


145.7 


7639 


2007-4-23 


ACIS-I 


00 h 25 m 22 s .02, 64°08'13" 


1 


108.9 



shock expected with no relativistic particles (e re i — 0%) 
are ~1.1 and ^1.3 times larger than those observed 
(e re i — 63%). This means that the observed m- value 
would be modified to 0.65 (= 0.54x1.3/1.1) when com- 
pared with the SNR evolutionary models without rela- 
tivistic particles. By using the m-value of 0.65, we derive 
the time to be 0.2 in characteristic units again based on 
Fig. 1 in Dwarkadas (2000). Then, the ambient density 

is calculated to be -O.OOlS^^Mcj/Mch) 572 cm" 3 . 
Since this model requires even lower ambient density 
than that from the model without particle accelera- 
tion (^0.2 cm -3 ), we conclude that the ambient density 
around Tycho's SNR is likely less than ^0.2 cm~ 3 . Fur- 
ther detailed SNR evolutionary models considering the 
effect of relativistic particles would be helpful in deriv- 
ing an accurate ambient density. 

5. CONCLUSION 

We have presented X-ray proper motion measurements 
of the edge (i.e., forward shock and ejecta knots) and the 



reverse-shocked ejecta of Tycho's SNR. The azimuthal 
variation of the proper motions along the edge as well 
as the average expansion index of ^0.5 are consistent 
with those derived from the most recent radio mea- 
surements (Reynoso et al. 1997). The proper motions 
for the reverse-shocked ejecta is measured to be 0".21- 
0".31yr _1 . Comparison of the expansion index of the 
forward shock and the reverse-shocked ejecta with an 
evolutionary model of a Type la SN gave us a typical 
pre-shock ambient density of less than ~0.2cm~ 3 . 
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TABLE 2 

Summary of proper-motion measurements for forward shock (FS) and ejecta knots (EK) 

Azimuth (degree) Proper motion (arcsecyr - 1 ) Expansion rate (%yr _1 ) Expansion index 

2000-2007 — 2003-2007 2000-2007 — 2003-2007 2000-2007 — 2003-2007 

3 (FS) 0.338 ± 0.016 — 0.256 ± 0.024 0.141 ± 0.007 — 0.107 ± 0.010 0.605 ± 0.029 

9 (FS) 0.374 ± 0.016 — 0.279 ± 0.021 0.160 ± 0.007 — 0.119 ± 0.009 0.687 ±0.029 

15 (FS) 0.309 ±0.011 — 0.270 ±0.015 0.135 ± 0.005 — 0.119 ± 0.007 0.582 ± 0.020 

21 (FS) 0.340 ± 0.015 — 0.273 ± 0.019 0.149 ± 0.007 — 0.120 ± 0.008 0.641 ± 0.028 

27 (FS) 0.325 ±0.017 — 0.261 ±0.024 0.142 ± 0.008 — 0.114 ± 0.011 0.612 ± 0.033 

40 (FS) 0.361 ±0.022 — 0.217 ± 0.029 0.156 ± 0.010 — 0.094 ± 0.013 0.672 ± 0.041 

45 (EK) 0.217 ± 0.019 — 0.239 ±0.021 0.085 ± 0.007 — 0.094 ± 0.008 0.366 ± 0.032 

56.25 (EK) 0.195 ± 0.016 — 0.213 ± 0.025 0.076 ± 0.006 — 0.084 ± 0.010 0.328 ± 0.027 

70 (FS) 0.263 ± 0.018 — 0.142 ± 0.029 0.117 ± 0.008 — 0.063 ± 0.013 0.503 ± 0.034 

78 (FS) 0.289 ± 0.013 — 0.173 ± 0.012 0.103 ± 0.006 — 0.079 ± 0.006 0.442 ± 0.025 

85 (FS) 0.307 ±0.011 — 0.204 ±0.012 0.108 ± 0.005 — 0.093 ± 0.005 0.466 ± 0.021 

92 (FS) 0.242 ±0.020 — 0.137 ±0.022 0.106 ± 0.009 — 0.060 ± 0.010 0.456 ± 0.039 

97.5 (EK) 0.219 ±0.010 — 0.231 ±0.013 0.086 ± 0.004 — 0.091 ± 0.005 0.370 ± 0.016 

107.5 (EK) 0.293 ±0.014 — 0.279 ±0.019 0.115 ± 0.005 — 0.109 ± 0.007 0.494 ± 0.023 

117.5 (EK) 0.311 ±0.012 — 0.302 ±0.019 0.122 ± 0.005 — 0.119 ± 0.007 0.525 ± 0.020 

127.5 (EK) 0.334 ± 0.015 — 0.328 ± 0.025 0.131 ± 0.006 — 0.129 ± 0.010 0.563 ± 0.025 

152 (FS) — 0.255 ± 0.046 

158 (FS) — 0.291 ± 0.031 

164 (FS) — 0.248 ± 0.030 

170 (FS) — 0.264 ±0.025 

176 (FS) — 0.267 ±0.022 

181 (FS) — 0.286 ± 0.025 

186 (FS) — 0.324 ±0.027 

195 (EK) — 0.364 ± 0.025 

206 (FS) — 0.363 ± 0.024 

212 (FS) — 0.344 ±0.023 

217 (FS) — 0.370 ±0.026 

223 (FS) — 0.355 ± 0.026 

229 (FS) — 0.369 ± 0.019 

235 (FS) 0.350 ± 0.014 — 0.380 ± 0.017 

246 (FS) 0.279 ± 0.013 — 0.354 ± 0.015 

251 (FS) 0.287 ± 0.011 — 0.368 ± 0.013 

257 (FS) 0.315 ± 0.010 — 0.377 ±0.013 

263 (FS) 0.393 ± 0.015 — 0.398 ± 0.020 

269 (FS) 0.315 ± 0.013 — 0.367 ±0.017 

274 (FS) 0.325 ± 0.010 — 0.366 ± 0.014 

280 (FS) 0.357 ± 0.019 — 0.362 ± 0.023 

286 (FS) 0.327 ± 0.013 — 0.359 ±0.019 

292 (FS) 0.328 ± 0.018 — 0.317 ±0.027 

298 (FS) 0.319 ± 0.012 — 0.361 ±0.017 

304 (FS) 0.276 ± 0.013 — 0.322 ± 0.015 

310 (FS) 0.314 ± 0.008 — 0.330 ±0.012 

317 (FS) 0.297 ± 0.015 — 0.254 ±0.021 

325 (FS) 0.254 ± 0.011 — 0.261 ±0.018 

337 (FS) 0.272 ± 0.022 — 0.249 ± 0.031 

Mean 0.302 (0.059) 0.121 (0.022) 0.522 (0.094) 

NOTE. — Errors quoted arc statistical 1-sigma uncertainties. Systematic uncertainties for proper motions, expansion rates, and expansion indices 
arc O^.Oeiyr -1 (or 0".072yr _1 ), ~0.024%yr -1 (or ~0. 028 % yr _ 1 ) , and ~0.103 (or ~0.122) for 2000-2007 (or 2003-2007), respectively. The 
values in blackets with mean values are standard deviations. 






127 


± 





005 





103 


± 





005 





106 


± 





004 





118 


± 





004 





149 


± 





006 





119 


± 





005 





126 


± 





004 





138 


± 





007 





127 


± 





005 





130 


± 





007 





130 


± 





005 





115 


± 





005 





134 


± 





004 





132 


± 





007 





113 


± 





005 





123 


± 





010 






107 


± 





010 





119 


± 





009 





119 


± 





007 





120 


± 





008 





114 


± 





011 





094 


± 





013 





094 


± 





008 





084 


± 





010 





063 


± 





013 





079 


± 





006 





093 


± 





005 





060 


± 





010 





091 


± 





005 





109 


± 





007 





119 


± 





007 





129 


± 





010 





101 


± 





018 





116 


± 





012 





098 


± 





012 





102 


± 





010 





101 


± 





008 





107 


± 





009 





120 


± 





010 





143 


± 





010 





135 


± 





009 





127 


± 





009 





134 


± 





009 





129 


± 





009 





134 


± 





007 





138 


± 





006 





131 


± 





005 





136 


± 





005 





141 


± 





005 





151 


± 





008 





139 


± 





006 





142 


± 





005 





140 


± 





009 





139 


± 





007 





126 


± 





011 





147 


± 





007 





134 


± 





006 





141 


± 





005 





113 


± 





009 





116 


± 





008 





112 


± 





014 






546 


± 





022 





444 


± 





021 





458 


± 





017 





507 


± 





017 





640 


± 





024 





513 


± 





020 





542 


± 





016 





559 


± 





031 





544 


± 





022 





559 


± 





031 





558 


± 





021 





494 


± 





023 





578 


± 





015 





567 


± 





029 





485 


± 





021 





528 


± 





043 






461 


± 





043 





516 


± 





040 





512 


± 





028 





517 


± 





036 





495 


± 





046 





407 


± 





055 





405 


± 





036 





361 


± 





042 





272 


± 





055 





341 


± 





024 





402 


± 





024 





260 


± 





042 





391 


± 





022 





473 


± 





032 





512 


± 





032 





555 


± 





042 





437 


± 





079 





499 


± 





053 





424 


± 





051 





441 


± 





042 





436 


± 





036 





463 


± 





041 





518 


± 





043 





617 


± 





042 





581 


± 





038 





551 


± 





037 





579 


± 





040 





555 


± 





041 





577 


± 





029 





138 


± 





006 





131 


± 





005 





136 


± 





005 





141 


± 





005 





151 


± 





008 





139 


± 





006 





142 


± 





005 





140 


± 





009 





139 


± 





007 





126 


± 





011 





147 


± 





007 





134 


± 





006 





141 


± 





005 





113 


± 





009 





116 


± 





008 





112 


± 





014 



TABLE 3 

Summary of proper-motion measurements for reverse-shocked ejecta 



Sector 


Proper motion (arcsec yr 1 ) 


Expansion rate (% yr 1 ) 


Expansion index 




2000-2007 — 2003-2007 


2000-2007 — 2003-2007 


2000-2007 — 2003-2007 


NE(0"-60°) 

SE(132°-180°) 

SW(180°-240°) 

W(240°-300°) 

NW(300°-360°) 


0.210 ± 0.005 — 0.222 ± 0.007 
— 0.266 ±0.013 


0.100 ± 0.002 — 0.106 ± 0.003 
— 0.133 ± 0.007 


0.429 ± 0.009 — 0.457 ± 0.014 
— 0.574 ± 0.029 


— 0.297 ± 0.008 


— 0.145 ± 0.004 


— 0.625 ±0.017 


0.291 ± 0.004 — 0.313 ± 0.006 
0.221 ± 0.004 — 0.237 ± 0.007 


0.138 ± 0.002 — 0.149 ± 0.003 
0.106 ± 0.003 — 0.118 ± 0.003 


0.596 ± 0.008 — 0.644 ± 0.013 
0.473 ± 0.008 — 0.509 ± 0.014 


Mean 


0.294 (0.040) 


0.143 (0.019) 


0.614 (0.084) 


NOTE. — Same as tabl 
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Fig. 1. — Chandra three-color image of Tycho's SNR after vignetting effects are corrected. Red, green, and blue represent 0.3-1.0 keV, 
1.0-4.0keV, and 4.0-8.0keV X-rays, respectively. The image is binned by 0".492 and has been smoothed by a Gaussian kernel of a = 0".984. 
The intensity scale is square root. The small boxes, small annular pieces with triangular marks, and a large annulus divided by five sectors 
are where we investigate the proper motion of the forward shocks, ejecta knots, and reverse-shocked ejecta, respectively. We show an 
example one-dimensional X-ray profile from the yellow box indicated by an arrow in Fig. [3] Example x 2 or C values calculated for the 
yellow box and the yellow annular sector with an arrow are shown in Fig. [4] and Fig. [6] 




Fig. 2. — Left: Difference image between 2000 and 2007 observations in 1.0-8.0 keV band. The original images before subtraction are 
binned by 0".492 and smoothed using a Gaussian kernel of a = 0".492. The intensity is linearly scaled from -3 to +3 counts pixel -1 . The 
southern dark (negative) region of the remnant is due to the fact that the observation in 2000 does not cover there. Right: Same as left but 
for 2003-2007. The structures running from NW to SE and from NE to SW are caused by the gaps of ACIS chips. Also visible artifacts 
are stripes running from NE to SW caused by bad columns. 
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Fig. 3. — Example one-dimensional profiles extracted from the yellow box shown in Fig.[T] Data points with black circles, red boxes, and 
blue triangles represent the 2000, 2003, and 2007 epochs, respectively. 




Shift (arcsec) 

Fig. 4. — Example \ 2 ~ an d C-value distributions as a function of shift position. The x 2 values with a solid line are calculated from the 
2003 and 2007 one-dimensional profiles shown in Fig.[3](see, section 3.1.1), while the C values with a dashed line are calculated from 2003 
and 2007 polar-coordinate images for the same region (see, section 3.1.2). 




Azimuth angle (degree, counterclockwise from N) Azimuth angle (degree, counterclockwise from N) 



Fig. 5. — Left: Proper motion for the forward shocks and the ejecta knots as a function of azimuthal angle (counterclockwise from 
the N). Black and red are derived from 2000-2007 and 2003-2007 measurements, respectively. Data points with filled boxes and open 
circles represent the forward shocks and the ejecta knots, respectively. The data for 2003-2007 measurements are displaced by 1 degree in 
azimuthal angle for clarity. The southern portion of the remnant (152°— 229° in azimuth angle) was not covered by the 2000 observation. 
Right: Same as left but for expansion rate (left axis) and expansion index (right axis). 
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Fig. 6. — Example C-value distribution as a function of shift position for the reverse-shocked ejecta from the yellow (NE) annular sector 
in Fig. [T] calculated from 2003 and 2007 polar coordinated images. 
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7. — Same as Fig. [5] but for the reverse-shocked ejecta. 
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